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We find the emergence of strong correlations and universality on the approach to the quantum crit- 
ical points of a two impurity Anderson model. The two impurities are coupled by an inter-impurity 
exchange interaction J and direct interaction U12 and are hybridized with separate conduction 
channels. The low energy behavior is described in terms of renormalized parameters, which can be 
deduced from numerical renormalization group (NRG) calculations. We show that on the approach 
to the transitions to a local singlet and a local charged ordered state, the quasiparticle weight factor 
z — > 0, and the renormalized parameters can be expressed in terms of a single energy scale T*. 
The values of the renormalized interaction parameters in terms of T* can be predicted from the 
condition of continuity of the spin and charge susceptibilities, and correspond to strong correlation 
as they are greater than or equal to the effective band width. These predictions are confirmed by 
the NRG calculations, including the case when the onsite interaction (7 = 0. 

PACS numbers: 71.10.Ay,71.10Hf,71.27.+a,73.63.Kv 



There is increasing interest, both experimentally and 
theoretically, in strongly correlated electron systems 
which have anomalous behavior in the region of a T = 
or quantum phase transition (for a recent review see [l[). 
The generalizations of the Wilson renormalization group 
approach to zero temperature transitions, initiated by 
the early work of Hertz and followed up by others [2j, 
have not provided a comprehensive framework to explain 
many examples of quantum critical behavior. Particu- 
larly challenging is the range of anomalies observed at 
quantum critical points (QCP) in heavy fermion materi- 
als, which have been induced by lowering the transition 
temperature of a magnetically ordered state to zero by 
pressure, alloying or in some cases by an applied mag- 
netic field P, 0, 0| • 

One possible mechanism that has been put forward for 
some heavy fermion systems is that, at the critical point, 
there is a breakdown of the Kondo screening such that the 
associated Kondo resonance disappears [l|,|j|. This would 
imply that the f-like quasiparticles no longer contribute 
the Fermi surface, so that the Fermi surface would shrink 
from a large to a small one, containing only the itiner- 
ant non-f electrons. It is difficult to test this conjecture 
for a lattice model of this situation using the theoretical 
techniques currently available. There is, however, a two 
impurity Kondo model which has been shown to have 
a quantum critical point . which might throw some 
further light on this particular mechanism. In this Letter, 
we study a related two impurity model using a combina- 
tion of renormalized perturbation theory (RPT) and nu- 
merical renormalization group (NRG) calculations. Us- 
ing these techniques we predict that a single energy scale 
T* emerges as the critical point is approached, such that 



at the critical point T* = 0, and the quasiparticle Kondo 
resonance disappears. These results lead to a new per- 
spective on the two impurity model, and evidence sup- 
porting the conjectures that a quantum critical point can 
be associated with the collapse of the Kondo quasiparti- 
cle resonance. 

The Hamiltonian of the model we will study has the 
form, % = X)a=i 2 ^« + ^12 where H a corresponds to 
an Anderson impurity model in channel a given by 

Ma = e d,ad\ y , rJ d a (T + € k,ac\ a iJ C k a (J (1) 

+ y^(Vfc,a4 t ,o- c fc,a.<T + n - c + U a n d ^ a ^n dta ^ 

where d) a CT , d a a , are creation and annihilation operators 
for an electron at the impurity site in channel a, where 
a = 1, 2, and spin component a =t, ^. The creation and 
annihilation operators ct , c k are for partial wave 
conduction electrons with energy ek, a in channel a. 

The second part of the Hamiltonian H12 decribes the 
interaction between the impurities in the two channels, 
which we take in the form of a direct Coulomb term U12 
and an Heisenberg exchange term, 

Hl2 = U12 nd,l,<T J! n d,2,a' + 2 JS d: i • Sd,2, (2) 
a a' 

where J > for an antiferromagnetic coupling. 

For the model with U\2 = there is a competition 
between two modes of screening of the impurity spins; 
by Kondo screening in the channel directly hybridized 
to each impurity or by the direct antiferromagnetic cou- 
pling between the impurities. The Kondo screening by 
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the individual channels predominates for small J, and 
the local screening for large J. NRG studies [j| of the 
Kondo version of the symmetric model have shown that 
there is a critical point between these competing terms 
at a value J — J Cl where J c is proportional to the Kondo 
temperature 7k of an isolated impurity (J = 0). It was 
also shown there is a divergence of the impurity specific 
heat coefficient at the transition point and an anoma- 
lous log(2) /2 entropy. Conformal field theory 0, [l(| and 
bosonization studies Q have clarified the conditions for 
such a transition to occur and shown that the staggered 
susceptibility diverges at the transition point. 

The model with J = 0, U\2 ^ 0, has been used to de- 
scribe two capacitively coupled quantum dots, and NRG 
studies have revealed a second type of transition at a crit- 
ical value of U\2 = U{ 2l such that for Un > U± 2 there 
is a breaking of local charge order Further recent 

NRG studies of both types of transitions have been re- 
ported using related versions of this model to describe 
the electron transport in double quantum dots Il2l-ll6l|. 

In the renormalized perturbation theory [17'), the 
impurity retarded Green's function Gd. a ,o{u) is re- 
expressed as G d , a ,<j{u) = z a G dta ^(uj). where G d , a ,a(u) 
is the quasiparticle Green's function given by 



Gd,a,a{u) — 
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l d + iA - Scr(w) 



(3) 



Exact results for the spin susceptibility and charge sus- 
ceptibilities [18] . Xs and Xa are given by 

Xs =Vl^(°)(0), f ]s = l + (U-J)p^(0), (6) 
Xc=47) cy 5(°)(0), fj c = l-(U + 2U 12 )p^(0), (7) 

and the phase shift S = tan _1 (A/erf) per spin per chan- 
nel. 

We consider first of all the symmetric model with 
U12 = 0. When J = the model corresponds to two 
independent Anderson models and the low energy be- 
havior corresponds to a local Fermi liquid in terms of the 
two renormalized parameters A and U. 

When we switch on and increase an antiferromagnetic 
coupling J this quasiparticle Fermi liquid picture breaks 
down at a particular coupling J = J c , when the quasi- 
particle weight z — » implying A — > 0. This in turn 
implies that the specific heat coefficient 7 given by Eq. 

diverges at this point. However, as we approach this 
point we would not expect the local uniform spin and 
charge susceptibilities to diverge as these susceptibilities 
are suppressed by an antiferromagnetic coupling and the 
critical fluctuations occur in the staggered spin channel. 
From Eq. ([7]) this can only be avoided if in turn the co- 
efficients fj s — > and r/c — > at the same point, so the 
product with the singular part remains finite. Assuming 
U12 = these conditions imply 



and the renormalized parameters, and A are given by 



e d = z(ed + £ CT (0)), A = zA, 



(4) 



where z = 1/(1 — <9£ CT (w, 0)/doj) evaluated at lu = and 
A = 7r J2k I 1 2 <5 (efc ) . We have taken the two channels to 
be equivalent and have dropped the channel index. 

In working with the fully renormalized quasiparti- 
cles, it is appropriate to use the renormalized or ef- 
fective interactions between the quasiparticles which 
we identify with the renormalized local four vertices 
z 2 r"'^, a „ cr „, {uji, u>2, 0J3, W4) in the zero frequency limit, 

eg. u a = z 2 r;;^ t (o, 0,0,0) 0. 

The effective Hamiltonian which describes the low en- 
ergy excitations corresponds to the original model given 
in equations ((I) and © with the parameters replaced by 
the renormalized values, and the interaction terms have 
to be normal ordered. For the complete renormalized 
perturbation expansion to include higher energy scales, 
counter terms have also to be included to cancel off any 
further renormalizations (l7| . 

The quasiparticle interaction terms do not contribute 
to the linear specific heat coefficient 7 of the impurities, 
which is given by 



7-47r 2 p (0) (0)/3, p {0 \uj) 



A/ir 



(5) 



(w - id) 2 + A2 
where p(°\ui) is the free quasiparticle density of states. 



J ->• 27rA, U -> 7rA, as J — > J c 



(8) 



We can define an energy scale T* via nA(U, J) = 4T*, 
such that evolves continuously from the value Tk for 
J = 0. This implies that on approach to the critical 
point the renormalized parameters can all be expressed in 
terms of this single energy scale, J/2 = U = ttA = 4TJj. 

We can apply precisely the same arguments to the 
model with finite Ui 2 and J = 0, which has a local charge 
ordered transition as Ui2(> 0) is increased to a critical 
value Ui2- Assuming J = 0, we find 



U 12 -> ttA, U -ttA, 



U12 -> U{ 2 



(9) 



Using the RPT approach [18j we can calculate the exact 
asymptotic behavior of the impurity retarded self-energy 
T) for ui, T « T< from the second order calculation 
of E(w,T) as ImE(w,T) = ImS(w,T)/z. The result is 



£(w,T) 



-i7T 2 7A 

"~64 



Ti 



ttT 



(10) 



where I = (2U 2 + 3J 2 + 4Uf 2 )/(nA) 2 , so I -» 14 as 
J ^ J c and I -> 6 as U -> C/f 2 • 

We have shown in earlier work how the renormalized 
parameters can be deduced from an analysis of the low 
energy fixed point of an NRG calculation [3 EH' and 
we apply the same procedure for this model to calculate 
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FIG. 1. (Color online) A plot of A/A, U/nA (coincident) 
and J/nA as a function of J/J c for U/nA = 5, nA — 0.01. 

the low energy behavior and test the predictions given in 
Eqns. ([5]) and In Fig. Q]we give results for the ratio 
of z = A/A. ?7/7rA and J/irA for £//7rA = 5 as a func- 
tion of J/J c , where J c = 1.378Tk- It can be seen that all 
three renormalized parameters tend to zero at the tran- 
sition point. The fact that z — > implies that the Kondo 
resonance at the Fermi level disappears at the transition. 
The ratio of parameters, U /ttA and J/ttA remain finite. 
For U 1 7rA = 5 we are in the strong coupling limit of 
the Anderson model (J = 0) where U/ttA — 1, and it 
remains at this value over the whole J interval. This is 
not the case in general, however, as can be seen in Fig. 
[2] (left) where values of U/nA are given as a function of 
J/J c for U/ttA — 0,0.5,1,2,6. It can be seen that as 
J —> J c , U/ttA — > 1 in all cases including the case U = 0, 
so that strong correlation result emerges even in the weak 
coupling case on the approach to the critical point. In 
the right panel of Fig. [5J the corresponding values of 
J /tt A are shown and all converge to the limiting value 
J /ttA — 2, which confirms the predictions given in Eqn. 
([5]). The results also apply for the model with U12 7^ 
as it is found that Uyi/'kA — > as J — > J c . 




FIG. 2. (Color online) A plot of U/irA (left) and J/irA 
(right), as a function of J/J c for U/nA — 0,0.5, 1,2,6. 

At J = J c there is a discontinuous change in the NRG 
fixed point as the two impurities decouple from the con- 
duction band on the lowest energy scale. The phase shift 
6 changes from tt/2 to due to the singularity develop- 
ing in the self-energy of the impurity Green's functions. 
For J > J c , z — and the previous analysis based on 
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FIG. 3. (Color online) A plot of the renormalized parameters, 
r</Tk, I>/Tk (stars) J/U (crosses), as a function of J/J c 
for U/tvA = 5. 




FIG. 4. (Color online) A plot of the Wilson ratio fj s (left) and 
the spin susceptibility ratio, Xs(U, J)/xs(U, 0), as a function 
of J/J c for U/ttA = 5. 



the assumption of analyticity of the self-energy at lo = 
breaks down. However, in this regime the low energy be- 
havior still corresponds to a Fermi liquid. We can retain 
Eqns ([7]) and ([5]) as a description of a local Fermi liquid 
and treat the first conduction site in the NRG chain as 
an effective impurity. We can then derive effective renor- 
malized parameters as for J < J c , but we have to take 
into account that the hybridization is now to a modified 
conduction chain. The renormalized quantity A is no- 
longer equal to zA, so in using it to define an energy 
scale T* we distinguish it from the values for J < J c , by 
AT* = irA. 

In Fig. [3] we show the results for the renormalized 
parameters over the range through the transition point 
as a function of J/J c for U/ttA = 5. It can be seen that 
J/U is continuous through the transition and takes the 
predicted value 2 at J = J c . The curves for T< and T> 
approach the critical point in a similar way proportional 
to ( J— J c ) 2 so that it seems reasonable to identify them as 
a single energy scale T* . The results shown were found to 
be universal in the strong correlation regime U /ixA > 3. 

The value of fj s , which is the Wilson ratio, is shown as a 
function of J/J c in Fig. 0](left). The corresponding spin 
susceptibility Xs = 4/ig?7 s /7rA, is difficult to determine 
precisely in the immediate region of the critical point as 
both A and fj s tend to zero in this limit. The value of 
fj s depends on the difference between the renormalized 
parameters which are very small in this regime, so any 
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errors in the determination of the parameters become 
significant. As J is increased to J/J c = 0.95 there is an 
almost linear decrease with J/J c . For J/J c > 1.05 the 
susceptibility appears to fall off more slowly with increase 
of J/J c . An interpolation between the two regimes, as 
shown in Fig. U (right), suggests that there could be a 
peak at the critical point J = J c , but it is very sensitive 
to the range of the interpolation regime. 
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FIG. 5. (Color online) A plot of A/A (left) and U/ttA (cir- 
cles) and Ui2/nA (stars) (right) as a function of U12/U in the 
approach to the charge order transition for U /nA = 5. 



Results for the transition to the local charge ordered 
state for the model with U12 ^ and J = are shown in 
Fig. [5] for U/irA = 5 as a function of U12/U, As U12 is 
increased there is a point U12 = U with SU(4) symmetry 
where U12 = U = 7rA/3, which is predicted from the fact 
that for large U both the spin and channel fluctuations 
are suppressed. The critical point occurs for U > U12, 
where A — > implying z — > and the disappearance 
of the resonance at the Fermi level. There is a rapid 
reduction in U/ttA from the SU(4) point to a value - 
1 at the transition and a commensurate increase in the 
value of l7 12 /7tA to the value 1, in complete agreement 
with the predictions based on Eqn. ([7]). As the quantum 
critical point (QCP) is approached we again have a single 
energy scale T* such that U12 = -U = AT*/tt, J = 0, 
which are also found to apply for the model with finite J 
as J/ttA — > at the transition. 

In summary, we see that universality appears on the 
approach to the quantum critical points, such that the 
renormalized parameters specifying the low energy be- 
havior can be expressed in terms of a single energy scale 
T*. At the critical points T* — > the quasiparticle weight 
factor z — > and the spectral density of the impurity 
levels at the Fermi level goes to zero. The quasiparti- 
cle interactions are equal or greater than the renormal- 
ized effective band width 7rA, as in the strong correlation 
regime. The arguments used here should be generally 
applicable to models of heavy fermions as all the suscep- 
tibilities in Fermi liquid theory at T — take the form 
Xa oc p(0)fj a where p(0) is the density of states of the 
non-interacting quasiparticles at the Fermi level and fj a 
is a factor which depends on the interactions between 
the quasiparticles. If the specific heat coefficient, which 
is proportional to ,5(0) diverges at the QCP and the sus- 



ceptibility Xa is finite, then fj a = gives a constraint on 
the quasiparticle interactions. The emergence of a single 
low energy scale T* means that the low energy dynamic 
response functions would have the form F(u)/T*,T/T*). 
This would be a natural precursor of w, T scaling because 
as T* — >• 0, it would be expected to go over to a form 
T 7 /(w/T, 1). Calculation of the renormalized parame- 
ters from the NRG within a dynamical mean field theory 
for a lattice model would require the self-consistent solu- 
tion of the effective band conduction density of states for 
a two band model. 

Some recent interesting experiments have set out to ex- 
amine the QCP in a two impurity Kondo model by mea- 
suring the current between a cobalt atom on an STM tip 
and a cobalt atom on a metal surface [20|. The results 
are given as a function of the bias voltage so are under 
non-equilibrium conditions. There is a direct hybridiza- 
tion term between the cobalt atoms which is not in our 
model but could be included. Once the renormalized pa- 
rameters have been determined it is possible to calculate 
precisely the differential conductance at low bias voltage, 
using the Keldysh version of the renormalized perturba- 
tion theory. This approach could be used to calculate 
the onset of the splitting of the Kondo resonance seen in 
these experiments in a similar way to the calculation of 
the onset of the splitting in a magnetic field in a quantum 
dot 
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